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Abstract 

This paper deals with global existence and uniqueness for the lake 
equations with a bottom topography vanishing on the shore. Our 
result generalizes previous studies that assumed the depth to be non- 
degenerate. Elliptic estimates for degenerate equations are established 
studying the behavior of the associated Green function. 
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1 Introduction 

In this paper we are interested in a two-dimensional geophysical model that 
has been essentially described by Greenspan in |Gr] page 235. This sys- 
tem describes the evolution of the vertically averaged horizontal component 
v{t,x,y) of a three dimensional velocity vector u{t,x,y, z) and reads 

dt{bv) + divx(6-u ^ v) + bV„p = in 17, 

div^r(bv) = in rj, 
(1.1) { , 

[bv) • n = on dQ, 

ibv)\t=o = mo in Q 
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where x denotes the horizontal components x = {x,y). The lake equations 
may be seen as the low Froude number limit of the usual inviscid shallow 
water equations when the initial height converges to a nonconstant function 
depending on the space variable, namely 6(x), see for instance [BrGiLij . 
Recall that the inviscid shallow water reads 



(1.2) 



' dth + divx(/if ) = in 17, 

dt(hv) + divjhv v) + Wx^^W^ = in 17, 

Fr^ 

(hv) • n = on dil., 
^ {hv)\t=o = mo, h\t=o = homn, 



where Fr is the Froude number, h is total depth, v the vertical mean value 
of the horizontal components of the velocity. 

The constant case b = const.- The 2D incompressible Euler equations. We 
remark that the case 6 ~ 1 reduces to the standart 2D incompressible Euler 
system. The Cauchy problem for incompressible Euler equations is very well 
understood and we refer the reader to various existing surveys on the ques- 
tion : see for instance [Yoj . |Ma] . |Ch| and |BMj . Concerning the bounded 
domain case, we refer to [Li] where the method used by Youdovitch is ex- 
plained and various results described. The approach is the construction of 
the solutions as the inviscid limit of solutions to a Navier-Stokes system 
with artificial viscosity and boundary conditions v ■ n = and uj = where 
u! = dyVi — dxV2 is the vorticity. We also mentionned |CMRj where a wall 
law related to the cauchy stress tensor replaces the null condition on the 
vorticity. In all these works, to get the result, elliptic estimates are used on 
the following equation 

-Ax^- = tj in 17, = 0. 

Namely the following estimate 

< ||w||ip(n) 

where C does not depend on p. Such elliptic estimates with nondegenerate 
coefficients follow from [AgDoNiT - |AgDoNi2 . We also mention in the 



method based on approximate solutions which strongly uses the following 
property on VxW, see for instance |BM] : Let vhea. smooth, L^(R^)nL°°(R^), 
divergence free field and let uj = curlu, then 

< c(l + ln+ |b||L3(]R2) + ln+ ||a;||i2™2))(l + \\u:\\l 
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where ln~^ g = ii g < 1 and In^ gf = In (7 if (7 > 1 . 

The nonconstant case b 7^ const; The lake equations. If the depth is not 
assumed to be constant, we have to define a new vorticity variable to get a 
transport equation on this quantity. Namely defining 

we get that it is transported along the flow v namely 

(1.3) dtuj + v VxW = 0. 

Together with the incompressibility condition div (bv) = 0, this implies, 
formally, conservation of the U'{bdx) norms of co, for p € [1, 00]. 

The case 6 > C > has been studied in [LOT] proving global existence 
and uniqueness for an initial vorticity lvq in L°°(f2) following Youdovitch's 
method. The regularity of the lake equations has been also addressed in 
[Ca] in this case but where |V6(3;)| ^ 00 as x tends to the boundary. Here 
we study the degenerate case, namely the case when b stricly positive in Q 
vanishing on the shore dQ. Note that, even if the model seems to be not 
physically relevant in the presence of beaches, similar models are used by 
several geophysicits. In a forthcoming paper, the authors will try to better 
understood the right model to be considered in the presence of beaches. 

We will assume that b > on Q and 6 = on dQ. More precisely, we 
assume that is bounded with smooth boundary d^l, and 

(1.4) b = ip'' 

where a > and Q = {ip > 0} with (p E C°°{Tl) and V^ip 7^ on dfl. 

In order to get our result, an elliptic regularity result is required for a 
degenerate equation on the associated stream-function. This is the main 
part of the paper. It is strongly related to a careful study of the associated 
Green function given in |GSj - [GST| . 

The paper is organized as follows. The main results will be summarized 
in Section 2. In section 3 we will mention the existence and uniqueness result 
and the key result of the paper: Namely the elliptic estimates on the stream 
function ^. Sections 4, 5, 6 are respectively devoted to the proof of such 
LP estimates by proving the normal regularity and vanishing property, by 
showing the Holder estimates and finally by establishing the estimates 
respectively. The well posedness result will be proved in the last section 
adapting Youdovitch's proof used for the standart 2D incompressible Euler 
equations. 

For the sake of simplicity, throughout the paper, we will supress the x 
index from the partial derivatives. 
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2 Main results. 



The goal of the paper is to prove an existence and uniqueness result of a 
global weak solution for the inviscid equations 11.11 We assume that is 
a smooth simply connected bounded domain and that b is given by ()1.4p . 
As usual, we eliminate the pressure p from the equation and consider the 
voriicity- stream formulation of the 2D lake equation. The weak form of the 
transport equation (jl.3|) is 

(2.1) dtibio) + dwibiov) = {). 

The vorticity and and the velocity are linked by the equations: 

(2.2) div ibv) = 0, curlu = boo, (bv) • n^gn = 0. 

This leads, see for instance [MaPuj . to introduce a stream function ^ such 
that 

(2.3) v = lv^^ = lid2^,-di^), ^\dn = 0. 

b b 

In the non-simply connected case, which is discussed later in a remark, one 
has to specify different boundary values on each connected component of 
dVt, see |MaPu| . The model for ^ then reads 

(2.4) div(iv^') = 6w, ^|an = 0. 

Note that, by Hardy's inequalities, the space Co°(J7) is dense in the space 
G Hl{^)]b-^/'^V^ G L'^iVt)]. Therefore, for / G ^^(Q) the problem 

(2.5) div(^V^) = /, ^|af7 = 0, 

has a unique solution ^ in this space, denoted by ^' = K f . We make the 
following definition. 

Definition 2.1. Given loq G L°°(r2), (f,u;) is a weak solution to the vorticity- 
stream formulation of the 2D lake equation with initial data wq, provided 

i) G L~([0,r]) X Q) andbuj G C^{[Q,T];L'^^{n)), 

ii) bv = V^K{bio) G C0([0,T];L2(n), 

iii) For all ip G C^{[0,T] x Q) and ti G [0,r].- 

/ bip{ti,x)uj{ti,x) dx— / bip{0,x)LL!o{x) dx = / / {bdt(p + bv ■V(p)u! dxdt. 
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Theorem 2.2. i) (Regularity) Assume that to G L°°((0,r) x Q) and 
bv G L°°([0, T]; L^(r2)) such that div{bv) = satisfy the weak formulation. 
Then cj G C°{[0,T]; U' (n)) and v G C°([0, T]; Tyi'''(0)) /or a// r < +00. 
Moreover, there is C such that for all p > 3, 

In addition, we get the following boundary condition on v 

V ■ n = on do,. 



ii) (Existence) For all ujq G L°°{0,), there exists a global weak solution 
{v,uj) to the vorticity-stream formulation of (jl.ip . 

iii) (Uniqueness) The weak vorticity-stream solution is unique. 

This result fohows the Youdovitch's procedure in constructing the so- 
lution as the inviscid limit of solutions of a system with artificial viscosity 
which is the analog of Navier-Stokes with respect to the Euler equations. 
The core of the proof is a regularity result for a system degenerate due to 
coefficient vanishing on the boundary of the domain. More precisely, the 
main part of the paper will concern the following main result. Let be a 
smooth simply connected bounded domain in M" with n > 1 and b defined 
by ()1.4p . Consider the system: 



(2.6) 

for 

(2.7) 



div (bv) = in ri, {bv) ■ n\Q^ = 0, 
curlu = / in O 



bv G L^{^), 

feLP{n). 



Using the definition of b given in ()1.4p . we prove the following result on 
which the existence and uniqueness result is based. 

Theorem 2.3. If{v,f) satisfy ^M) and dMD with p £\2,oo[, then 

(2.8) V G Ci-"/P(n), Vv G LP{VL). 
There is a constant Cp independent of {u,uj) such that 

(2.9) \\v\\ci~^/p(Q) < Cp{\\f\\Lp + \\bv\\L2). 
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In addition 



(2.10) 



V ■ n = on di}. 



Moreover, for all pQ > 2, there is a constant C independent of {v,f) and 
p £ \pq,oo[ such that 



In this statement, for p g]0, 1[, C^{i}) is the usual space of continuous 
functions on which satisfy the Holder condition of order /x. In particular, 
([^ imphes that 



Remark 2.4. It is remarkable that the final W smoothness of v is inde- 
pendent of b in ()2.6p . When p = 2, it is still true that v G W^'"^ and indeed 
estimate ()2.1ip is uniform for p G [2,oo[. This is shown using estimates 
of |BGl iBCj in the proof below in place of those of [BCM] . However, in 
this paper we are mainly interested in these estimates as p — > oo, and we 
leave aside the case p = 2. Note also that (j2.9p follows from (j2.1ip and the 
Sobolev embedding W^'"^ C for = 1 — 2/p > 0. The proof of the uni- 
form estimate ()2.1ip given below relies on proving (j2.9p for a given po > n 
first, and on winning the final estimate for Vv, then. This is why we 
have stated them independently. 

Remark 2.5. We include here some remarks on the link between the 
vorticity-stream formulation and (jl.ip . 

1. Suppose first that v G C^{[0,T], H^(Q)) is a solution in the sense of 
distributions of (jl.ip . For u £ H^, the following identity holds: 



implying that co = cur\v satisfies transport equation (|2.ip in the sense of 
distributions thus condition Hi) in Definition [2Tj Moreover, (j2.2p is satisfied 
by (II. ip and by definition of lo. Thus v = 6~^V"'"'I', (see the beginning 
of section 3) and the stream function ^' G Hq satisfies ()2.4p . hence v = 
b-^V^K{huj) as in ii) of Definition EH 

2. Conversely, suppose that (w, v) is a weak solution of the vorticity- 
stream formulation. By part i) of Theorem 12.21 v G C'^([0, T], L^(r2)) and. 



(2.11) 




(2.12) 
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because w satisfies ()2.ip . we also have dtUJ G L~([0, T], In partic- 
ular, 

dw{bdtv) = 0, cml{dtv) = bdtuj G L^{[0,T], H-^{n)), 

implying that dtV G L°°{[0,T], Lf^^{Q)). The same computation as in 1), 
implies that in the sense of distributions 

curl (^dtv + ^ div {bv ^ v))^ = dt{buj) + div (bvuj) = 0. 

Therefore, assuming that Q is simply connected, there is a pressure p G 
L°°([0,r],P'(17)) such that is satisfied. 

Remark 2.6. When has some islands that means it is a non-simp ly 
connected domain, some generalized circulations are specified in order to 
uniquely determine the velocity v, see |MaPu| . More precisely let O with n 
islands Ji, . . . , In, the missing parameters are the boundary values Aj with 
i = 1, . . . , 71 of the stream function on each of the boundary component dli, 
. . . , din- Assuming the generalized circulations to be zero, the velocity is 
given hy V = vq + "^^^i with Vj = (V^^'j)/6 where '^j {j = 0, . . . , n) 
are the unique solution of 

-div(^V^'i) = 0, mi\Qj^=8i^j j = 0,...,n i = l,...,n. 

and 

-div(-V^'o) = bio, ^-obc = 0. 



The n coefficients Xj are uniquely determined using the zero generalized 
circulations conditions. This decomposition shows that our study concerning 
the simply connected domain may be applied to the non-simply connected 
one, the estimate remaining true. 

3 Localization 

Since bv G L^(r2) satisfies div (bv) = and (bv) • n = on dCl, there is a 
unique potential ^' such that 

(3.1) v = ^V^^ = ^{dy^,-d,^), ^|9o = 0. 

Indeed, ^' is determined in ffQ(r2) as the solution of 

(3.2) -A^'= curl(6t;), ^-19^ = 0. 
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There holds 

curl (bv) = b curl v + \7b x v = bf — ^Vb ■ 
With (jl.4p . we obtain the equation 

(3.3) - vjA^r + aV(^ • = = 0. 

At this stage, the dimension d = 2 plays no particular rule and we can 
consider this equation on any smooth bounded domain = {99 > 0} C M", 
with (p £ C°°(M") and dip ^ on dil.. One introduces the Holder spaces 
C^{Q) and the Sobolev spaces W^'''^{Q) based on L^, together with the 
weighted spaces Cf (O) and W^'P{n) of functions u in C^(0) or VF^'^(r2) 
such that (fu in C+^Q) or W''+^^P{n) respectively 

Writing 

= (^"+1$, u = ip-^'V-^^ = ipV^^ + (a + l)^V-^ip, 
Theorem 12.31 follows from 

Theorem 3.1. // ^' G H^{^) and f £ 1^(9) with p > n satisfy ([33]) . then 
$ = (^-("+1)^ satisfies 

(3.4) $ G fl (H), $ G fl Tyi^'P(O). 

tJ.<l-n/p q<p 

Moreover, there are constants C^, ^ G [0, 1 — n/p] such that for all ^ G 
and f G LP satisfying ()3.3p . 

(3.5) \mc^(n)<CM\\L^ + \\nm), 

For all pq > n, there is a constant Ci, such that for all \I' G Hq and f £ LP 
with p > po satisfying (|3.3p 

(3.6) ^ll^llH^i.P<t^l(ll/l|L. + ||^||Hi). 

Equation ()3.3p is a particular case of degenerate elliptic equations studied 
in \BG\ \BC\ \GS\ IBCMj and the estimates (|3.5p (13. 6p are mainly consequences 
of the results obtained in these papers. Because it may be useful in other 
circumstances and also because it helps to understand the analysis, we will 
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prove the estimates for a slightly more general class of equations, namely 
equations of the form 

(3.7) L{x,d^)^ :=-ip{x)P2{x,d,)^ + aPi{x,d,)^ = ^''+'f 
where a is a positive real constant and 

n n 

(3.8) -^2 = ^ Pj,k{x)djdk: Pi = ^Pj{x)dj 

j,k=l 3=1 

have real and smooth coefficients on with pj^k{x) = Pk,j{x). We further 
assume that P2 is uniformly elliptic on $7 

(3.9) 5^Pi,fc(2;)e*efc > for X G 17, e G M"\{0}, 
and 

(3.10) Pi(x,0 = ^Pj,fe(x)Cj(9fc99(x) ondn. 

Notice that (j3.9p and (|3.1U|) are independent of the choice of the defining 
function ip for 17. 

Theorem 3.2. With assumptions as above, if G Hq(Q) and f G LP{Q) 
satisfy (13. 7p . then <I> = (/p'^^+^^^f satisfies ()3.4p and there are constants 
such that the estimates ()3.5p and ()3.6p hold. 

Reduction to a neighborhood of the boundary. Standard elliptic theory 
implies that on any open subset ili dl 17, ^' belongs to VF'^'^(17i) for all p 
finite and thus to C''~^^(17i) for all ^ < 1 — n/p. Moreover, 

(3.11) m\c^+^{n,)<C^{\\f\\Lv + \\nm) 
and, for p G [2, oo[: 

(3.12) ^W^^niv^n,) < Ci{\\f\\Lv + wnm) 

Therefore, if x e C°°(n) is such that x - 1 G C'^(^^), then 

(3.13) L(x^) = v5"+^x/ + 9 ^|an = 0, 

where the commutator g is C^~^ by (|3.1ip and vanishes on a neighborhood of 
do.. Therefore, one can factor out any power of writing g = ^p"'^^ fi with 
/i G with norm bounded by the right hand side of (j3.11|) . Therefore, 
it is sufficient to prove the estimates (|3.5p - (|3.6p for functions <1> that are 
supported in an arbitrary small neighborhood of 917. 
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Local coordinates near the boundary. The boundary dO, is a closed smooth 
manifold. Consider a coordinate patch x' i— > j{x') from an open set oj C 
M""^ to dCl, with 7 G on tJ. Taking v{x') to be the inward normal to 
dQ at 7(x'), we parametrize a neighborhood of 7(w) by (x', £ ujx] — 5, 5[ 
by considering the mapping 

(3.14) T : {x',xn)^^{x') + xniy{x'). 

In these coordinates, the equation (j3.7p is transformed to 

(3.15) L^- = -x„P2* + aPi* = <+V- 

The ellipticity property (j3.9p is preserved as well as (j3.10p which now reads 

n 

(3.16) Pi{x',0,O=Y.Pi'n{x',0)Cr 

i=i 

Collecting all these remarks, we see that Theorem 13.21 follows from the 
next estimates. 

Theorem 3.3. If^ G H'^{ujx]0,6[) and f G LP{ljx [0,6]) with p > n satisfy 
(I3T3D and 

(3.17) ^-,,=0 = 0, ^|{x.>V2} = 0. 
then ^ = x^+^<I> with 

(3.18) $ G C7|""/^(cJi X [0, (^]), $ G 'P(wi x [0, 6]) 

for all relatively compact open subset loi C uj. Moreover, there are constants 
Cfj_, /u G [0, 1 — n/p\ such that for all ^ G Hq and f G LP , p > Po > n 
satisfying ()3.3p . 

(3.19) ll<f|lcr{^,x[0,5[) < C-Mdl/llL^ + ll^llHi), 

(3-20) Imwl-^i^^.m) < Ci(||/||^. + 11^11^0- 

Here Cf [resp. W^'^] denote the spaces for functions u G [resp. 
u G VF^'*'] such that x„m G (7^^+^ [resp. XnU G VF^'^']. 
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4 Smoothness in normal direction and vanishing 
at the boundary 

The first step in the proof of Theorem 13.31 is to factor out x'^^ in ^' and to 
show that $ is Holder continuous in the normal variable. For simplicity, we 
drop the tildes in equation (j3.13p . 

Proposition 4.1. With assumptions as in Theorem \ 3.3\ there are fi > 0, 
C and a Banach space E C D^u), such that the function $ = Xn^"'^^^'^ 
satisfies on uj x [0, 6] : 

(4.1) cI>gC^([0,5];i?), a;„9,„<DGC7^([0,(5];i?). 

together with the estimates 

Multiplying the equation p.l3p by an appropriate factor, one can always 
assume that the coefficient of is exactly x„. Then, (|3.16p implies that 
the coefficient of dx„ in Pi is 1 when Xn = 0. Thus (j3.13p reads 

(4.3) -Xndl^^ + adx„^ = x^+^f + F, F = L'^, 
with 

(4.4) L' := ^ XnPj,kdjdk + Xn'^Pj,ndjdx„ - a'^pjdy 

j,k<n j<n j<n 

where the pn is smooth. 

We consider functions of Xn valued in (Banach) spaces of distributions 
in x' : L'^{[0,6];E) , C''{[0,6],E) etc , with E C P'(w). The space E is 
subject to change from line to line, and we note L^(P(g) , C^(P(g) when the 
space is unspecified. We also denote by x"P(g the space of products of the 
function x" with arbitrary distributions in T)'{oj). 

We know that 

(4.5) ^ (^H^ d C^I'^{L^). 

Using the equation ()4.3p . we also see that Xnd'^J^ € L^(iJ^^) and hence 

(4.6) xndx^^ G H\H-^) C C^'\H-^). 

We will make use of the next (classical) results (See also Proposition 4.1.2 
in [BCM] for more general results). 
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Lemma 4.2. For a > and fx g]0, 1[, the operators 



(4.7) IaU{Xn) = 




map continuously C^([0,5];£;) to C^([0, 5]; £;). 

Proof of Proposition \4.1\ Because ^' = for x > 5/2, Equation (j4.3p implies 
that 

(4.8) d,^ = Ux^^+^f + F)=x^J + UF) 
with 

(4.9) ^ = " / ^ Ci/^(LP). 

Because ^ € vanishes at x„ = 0, there holds 

rx„ 

(4.10) ^= / 5,„^'(t)dt = x„Ii(5,„^'). 

Jo 

a) Suppose that for some a £ [0,a[ there is fi €]0, 1[ such that 

(4.11) ^ex^C>^{V[^), x^d,^^ e x'^C'^{V[^) 

From (|4.5|) (j4.6|) . this is true when q = and fi = 1/2. 

This condition implies that F = x'^Fi with Fi S C'^(P[g). Therefore, 
J„(F) = x°Ja-a(Fi) and (gSI) implies that 

S.^M' = x-e + x^J,_,(Fi) e x"Ci/2(P[g) + <C^(P(g) C x^C^K^O 

for some /Ui > 0. This shows that (|4.11|) holds with a replaced by a + 1, 
decreasing fj. if necessary. 

Moreover, if (|4.11|) holds for some a, it also holds for all a' < a, with 
possibly smaller //'s. Therefore, after a finite number of iterations, one 
obtains that (j4.1ip is satisfied for some a £]a, a + 1[ and /.i > 0. 

b ) The property (j4.1ip implies again that F = x"Fi with Fi G C''(P(g). 
We now write 

Ja{F) = - x^/«_„(Fi) + xl / t^~^-^F^{t)dt 

Jo 

ex':,C^^{V')+xlV'cixlC^^{V') 
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for some > 0. Here we have used that a > a. 

With (113D this imphes that d^^^^ = x^*i with G C^'^(V'^^). Plugging 
this equation in (j4.10p . imphes that 

(4.12) ^ = xl+^Ia+ii-^i) G <+iC'^H^(g)- 

Summing up, we have proved that <I> = Xn^"'^^^'^ G C^^{T>[^) and 
xA^^ = - (a + l)^ G C^''{V[^)■ 

c) Inspecting the proof shows that the index ^ and the Banach space 
E such that (14. ip holds are independent of ^ and /. Moreover, each step in 
the preceding proof can be converted into an estimate ending with (j4.2p . □ 

5 Holder estimates of 

We continue the proof of Theorem 13.31 proving the Holder smoothness of $ 
and the estimate (I3.19P of Theorem 13.31 

Proposition 5.1. Suppose that ^ and f € LP([0,(5] x uj), p > n, are as 
in Theorem \3.3[ Then for all relatively compact subset toi C uj and all 
fi<l- n/p, ^ = Xn^"^^^^' G CiuJi X [0,5]) and Xn^ G C+^uJi x [0,6]), 
with 

(5.1) ||^'||c''(ZJix[0,5]) + lkn^'||cM+i(cJix[0,5]) < c{\\f\\LP + ll^ll/fi)- 

The equation (j3.15p for ^ implies that ^ satisfies 

(5.2) := -XnP2{x, 9x)$ - Qi{x, 9,)$ + R{x)<^ = f, 
where 

Qi{x,0 = 2(a + l)^Pj,„(x)^j - a'^pj{x)^j, 

j j 

R{x) = a{a + l){pn{x) - Pn,n{x)) / Xn- 

From (j3.16p . it follows that R is smooth up to the boundary {x„ = 0} and 
that 

n 

(5.3) Qi{x,0,0 = (a + 2)^p,,„(x',0)e,. 
Without loss of generality, we can also assume that = 1. 
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Before proving the proposition above, we introduce several notations. 
For fj, > 0, fi ^ N, we denote by C^(R_,_) the space of bounded functions 
on R" which are uniformly Holder continuous with exponent fi. Next, we 
denote by C^(l+), k £ N, the space of functions u £ C^(M_^) such that 
XnU £ C^~^^ for all j < k. Cl^ is the intersection of the C^. We refer to 
[BCM] for a definition of these spaces involving a Littlewood-Paley analysis 
and their extension to negative /u and real k. We also refer to the same paper 
for a precise definition of the spaces C^''^(M_,_), where the second index z/ G M 
measures the additional tangential smoothness. Finally, we introduce the 
spaces C7f "~ = C^''' ■ 

Proposition 15.11 is a direct consequences of the results proved in [BCMj , 
applied to second order equations (|5.2p . There are three assumptions in 
[BCM] to satisfy. 

• (HI) The second order part P2 is strongly elliptic. 

• (H2) The indicial polynomial, XP2{x' ,0,0,1) + Qi{x' ,0,0,1) has no 
roots in the strip Hi < Re X < fi2- 

• (H3) The differential operator = XnP2{x' , 0, irj, dx„)+Qi{x' , 0, irj, dx„) 
is injective in the space C^^(M+) for all t] G M"~^\{0} and all x' . 

1) For equation (|5.2|) the ellipticity assumption (HI) follows from (|3.9|) . 

2) By ()5.3p . the indicial polynomial is (A + a + 2). An important feature 
is that its roots are independent of x' . Since a > 0, (H2) is satisfied for all 

-1 < /Xi < fl2. 

3) Using again (j5.3p . we have 

(5.4) L° = XniZ^ - p^) + (a + 2)Z = Z{xnZ) + (o + 1)Z - x^p^ 
with 

dxn J]]r?jPj,„(0,x'), 

j<n 

^ Pj,k{x',0)r]jr]k. 

j,k<n 

Since p > when rj ^ 0, the bounded functions u in the kernel of are 
smooth and rapidly decreasing at infinity. Moreover, they are smooth up 



Z = 

(5.5) 

p' = 
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to 0, by the classical analysis of Fuchsian singularities. Thus, the following 
integration by parts are jusitified: 

= 3f?/ L\udx = - x„(|Zn|2 + pV)(ixn ^H0\^. 

Jo Jo 2 

This implies that u = 0, hence that the spectral condition (H3) is satisfied 
if /i + 1 > 0. 

Therefore we are in position to apply the results of [BCM] . 

Proof of Proposition 15. li First, we note that, because p > n, the Sobolev 
embedding theorem implies that the right-hand side of ()5.2p satisfies 

(5.6) /eLf cC-"/p(l^x [0,5]). 

a) <I> vanishes for x„ > 5/2, so we can extend it by for x„ > S. By 
proposition 14.11 we know that there are /u' s]0, 1/2] and a Banach space 
E C V'{lo) such that ^> G Cf'(l+;^). Therefore, Proposition 2.3.2 of 
|BCM] . implies that 

if X £ C°°{ui) is equal to 1 on oji. This implies that there is a (large negative) 
integer v such that 

(5.7) x^^^C^'^^'^iK), 
if /ii = /u' - 1 > -1. 

b) Let xi £ C!q°{u;), equal to 1 on toi and such that x = 1 on a neighborhood 
of the support of xi- There holds 

(5.8) Cxi^ = xif + [C,Xi]x^ 

By (ILSD, Xif e C/^^H^") C C^'-^'°{Wl) with /i = 1 - n/p. By (IST]) . 
[£, xi]x^ e C''i+i'''(l" ). Thus the right hand side of dSS]) belongs to C^^^" 
with = A* — 1 < < /xi + 1 = Since the assumption (H2) is satisfied 
for the pair 

(/ii,^2), Theorem 4.1 in [B(Mj implies that xi^' e ' °°, 

that is 

(5.9) xi^GCp+'''^^(M+), 
for some z^i. 
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c) Let X2 £ C^{io), equal to 1 on uji and such that x = 1 on a neighborhood 
of the support of X2- There holds 

(5.10) Cx2^ = X2f+[C,X2]xi'^. 

By <m, Xif G C^'-HK) C B y dsn), G C'^2+1''^i(R" ) C 

C^2''^i+i(M^). By Theorem 5.2 in [BCM] . we obtain that X2^ G (7^2+1,^2^ 
with 1/2 = min(0, i^i + 1). 

Therefore, after finitely many iterations, we conclude that there is x* G 
Cq°{lj), equal to 1 on such that 

(5.11) x*^ G Cf +^'°(M:^). 

d) Finally, introducing Xb G C^(u;), equal to 1 on toi and such that X* = 1 
on a neighborhood of the support of Xb writing the equation for Xb*^) 
Theorem 4.2 in [BCMj implies that Xb^ G Cf +^(1+) = Cf (1+) implying 
the proposition. □ 



6 LP estimates 

In this section, we finish the proof of Theorem l3.3^ proving the estimates 
(j3.20p . We recall that regularity result is performed in the general setting 
Q C with n > 1. Given ui as in the theorem, consider lv2 relatively 
compact in oj such that ZJi C a'2- By Proposition 15.11 <I> is of class on 
UJ2 X [0,(5] for /i < 1 — n/p. Consider x G C^(w2) with x = 1 on loi. Let 
= G C'f 1+). Then 

(6.1) CcP = g:=xf+[C,x]^(^LP{Rl). 

Moreover, by Proposition 15.11 choosing pQ > n and /xq = 1 — n/po, there 
holds for / G L^, p > pq: 

(6.2) Mc-^^iWl) < C{\\f\\L.o + \mm), 

(6-3) \\9\\lp{ri) < C'dl/IU" + ll^llcfO(S:^))- 

Therefore, the U' estimates (j3.20p follow from the next result. 

Theorem 6.1. Let /i g]0, 1[. Suppose that (j) £ Cf M^) has compact support 
in LUX [0,(5[ and C(p = g e LP{RI). Then (f) G Wlf^^i^^). 

Moreover, there is a constant C such that for all such cf) and p G [2, cx3[ 

(6.4) \\xndjdk(t)\\LP{wi) + \\dj(t>\\LP{wi) - ^{pWdWir-iwi) + ^WcA- 
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6.1 Preliminary results 

In this section we recall some known results about C alder on- Zygmund op- 
erators ([St]). We consider operators T acting in LJS^^p(M^) (the space of 
bounded functions with compact support) with kernel r(x, y) locally inte- 
grable away from the diagonal {x = y}. 

Proposition 6.2. Suppose that the kernel K{x,y) satisfies on x M"; 

(6.5) \Kix,y)\ < \d,Kix,y)\ < 

\x — y\^ ^ \x — y|" 

Then, the operator 

Tf{x) = j K{x,y)f{y)dy 
acts from L^comp^^^) to C^(M" ) for all fi < 1 — n/p. 

Proof. The first estimate for K implies that T maps LJ^^p(M^) to L°°. 
Moreover, the second estimate implies that 

I X — x' I f 1 

\K{x,y) - K{x',y)\ <C- for \x - x'\ < -\x - y\. 

\x — y\"' 2 

Thus, interpolating with the first estimate yields 

(6.6) \K{x,y) - K{x',y)\ < C———^ for |a; - x'| < -|x - y|. 

\^ y\ ^ 



We write 



T/(x) - Tf{x') = [ {K{x, y) - K{x', y))f{y)dy 

J \x-y\>2\x-x'\ 

+ [ {K{x,y)-K{x',y))f{y)dy. 

J\x-y\<2\x-x'\ 

By ()6.6p the first integral is 0{\x — x'\^). Note that for |x — y| < 2\x — x'\, 
there also holds \x' — y\ < 2>\x — x'\. Therefore, the second integral is 

0( I ^lL\=o(\x-x'\). 

^J\y\<^x-x'\\yY'-^^ ^ ' 

□ 
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Proposition 6.3. Suppose that T is a bounded operator in L^(M") with 
kernel K{x, y) satisfying for x ^ y : 

(6.7) \K{x,y)\ < , ^ , ■ \d^K{x,y)\< ^ 



|X — y|" ■ " ■ ■ \x — y\"'~^^ 

Then T maps LP{WI) to LP(M^) with norm 0{p) for all p e [2,+oo[. 

Proof. The adjoint operator T* is bounded in and its kernel K*(x,y) = 
K{y, x) which therefore satisfies 

(6.8) \K*ix,y)\<—^, \dyK*{x,y)\< f 

\x y\ \x y\ 

The Calderon-Zygmund theory imphes that T* is bounded from to 
the space weak-L^ (see e.g. [St]). Therefore, by Marcinkievic's interpolation 
theorem, T* maps to for p G]1,2] with norm 0{l/{p — 1) (same 
reference). By duality, this implies that T maps to for p £ [2, oo[ with 
norm 0{p). □ 

6.2 Parametrices 

Recall that C « —XnP2 — Qi, see (|5.2p . For y £ U, we denote by Cy the 
operator 

(6.9) Cy{x,d^) := -XnP2{y,d^) - (a + 2){d^^ +^Pj,niy)^j) 

j<n 

As above, we have assumed as we may that pn^n = 1- For a given y, there 
is a linear transformation 

(6.10) X = T{y)x, with Xn = Xn 
such that, in these variables Cy is transformed to 

(6.11) £ = -x„Ai- (a + 2)9£„. 

According to [GS] . Lemme 1, the fundamental solution of C is 

(6.12) E{x,y) = [ F{i,y, 

Jo 

with 

Hi, y,o) = i (y„)°+i^-("+") {e{i - e)r" , 
A^{i,y,e) = eD'^ + {i-e)b'^. 
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with 7 some constant depending on a and 

= |x„ - yn\ +\x -y = |x„ + yn\ + \x - y p. 

More precisely, for e g]0, 1[, let 

(6.14) E%x,y) = C ' F{x,y,e)de. 

Jo 

This function is singular only on y„ = and there holds 

(6.15) L{x,d^)E'{x,y) = G%i,y), 
where 

(6.16) G'{x,y) = 2{n-2 + ah y^+^^-^^+^+n) ^^^-^ _ (a+2)/2_ 
According to [GSJ , Appendix A, there holds : 

Lemma 6.4. G^{x,y) is an approximation of the identity as e tends to zero: 
it is nonnegative, converges uniformly to on compacts o/R^ xR^\{x = y} 
and for all bounded open set C Wl and all x £ Q 



lini / G%x,y)dy = 1. 



From E, we derive fundamental solutions for the operators Cz- 

(6.17) Ez{x,y) = \detT'{z)\E{T{z)x,T{z)y) 
and their approximate versions 

(6.18) EUx,y) = \detT'{z)\E'{T{z)x,T{z)y) 
Finally, we define the parametrices 

(6.19) E{x, y) = Ey{x, y), E'{x, y) = E^ix, y). 
Similarly, we define the Gl{x,y) and G^{x^y). 

Proposition 6.5. Suppose that a > 1. For G N, there is a constant Ck 
such that for {x,y) G x uj, x ^ y, e g]0, 1[, there holds 

(6.20) \V',E%x,y)\<- 

\x — y\"^'^ 

and for k > 1 

(6.21) \^nV'^E^{x,y)\<j—^^^. 
The same estimates hold for E. 
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Proof. It is sufficient to prove the estimates for E^{x,y) and therefore for 
E{x, y). For simphcity, we drop the tildes in the proof below. On M" x , 
there holds 

(6.22) A>\x-y\, A > {1 - e)^/\xn + yn). 
Hence, for a > — 1, 

^a+l^-{a+n) < 1^ _ _ 0y{a+l)/2^ 

implying 

<7|x-y|^-"r/2(l-0)-V2. 

Thus, integrating in 9, 

(6.23) |^^(x,y)| < C7|x-y|i-" 

Estimates on the derivatives of E. Differentiating k times A^, one proves 
by induction that 

(6.24) \V^^A\ < CkA^-''. 
Thus 

(6.25) |V^F| < (^(1 - 9))"/^ . 
Using (16:22]) yields 

I Vf;F(x, y)\ < C\x - y|i-'^-"r/2(l - 9)-^/'^, 

implying 

\VlE^{x,y)\<C\x-y\'-^-\ 
Sharper estimates of XnV^E. There holds 

(6.26) A^ = \x - y\^ + 4{1 - 9)xnyn. 
Let us assume first that 

(6.27) xr^yn < 2\x - y\\ 
Then, 

(6.28) {xn + ynf < {xn - Vn? + 4x„y„ < C\x - y\^ . 
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In this case, we use that A > \x — y\, without loosing any information and 
Integrating in 9 yields 

I XnV'Mx ,y)\< CXnVl^^ I X - . 

With ([OH]) . x„<+i < C\x - therefore 

(6.29) \xrylE'{x,y)\ < C|x - 
Consider next the case 

(6.30) XnVn > 2\x - y\\ 
In this case, 

(6.31) 5x„y„ > 2x2 + 2y2 ^ L<^<c. 
Let 

,:=^elO,ll. 

We use the estimates 

^2 > |x — when 1 — ^ < p, 

^2 > 4^]^ _ O^XnUn when 1 - 6* > p. 

Therefore, 

|x„VjE'(x, „)! < C-^i^l^ /'^ (»(1 - S)r"d6 
Thus 

|x„V^£;^(x,y)| < C , """^"t,^ p"/' 

(x„y„)("+'=+-)/2^ 
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For the last estimate we have used that n + k > 2. Therefore, 



X - y\n+k+a (-2,^y^)(a+2)/2 

T (t 11 \(n+k-2)/2 

XnUn \Xnyn) 



+ c 



a/2 

(6.32) \xnVlE'{x,y)\<C 



1, ' 1 

k cne/ M ^ r-i Vn J- 



a/2 Ij. yln+fc— 2 ' 

Together with (I6.3ip this imphes 
(6.33) \xnVlE'{x,y)\<C ^ 



This finishes the proof of the proposition. □ 
Lemma 6.6. The following identity holds 

(6.34) C{x,d^)E'{x,y) = G%x,y) + K'{x,y) 
where the kernels satisfy uniformly in e g]0, 

(6.35) |K^(x,y)|<^-^, \d^K^{x,y)\<^^^. 

Moreover, the kernels converge uniformly on compacts of {x / y} to 
K{x, y) which satisfy (j6.35p outside the diagonal. 

Proof. The definition of E^ imphes that it is singular only on y„ = 0. More- 
over, 

C{x, d^)E'{x, y) = Cy{x, d^)E%x, y) + K'{x, y) = G'{x, y) + K'{x, y) 
where 

K^ix^y) = Xny^{pj,k{y) -Pj,k{x))9x,0^^E'' 

(6.36) ^ 

+ ((a + 2)p,- „(y) - qj(x))d^^E' + R{x)E'. 

where the qj are the coefficients of Qi. In the first sum, the coefficient of 
Xndxjdxf.E^ is 0{\x — y\). By (j5.3p . the coefficient of dx^E"^ in the second 
sum is 0{xn) + 0{\x — y\). Together with Proposition 16.51 this implies the 
estimates (j6.35p . 

From (j6.36p . it is clear that converges to a kernel K which satisfies 
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Lemma 6.7. The kernels G^{x, y) are non negative, converge toO as e tends 
to uniformly on compacts subsets of (M_^ x uJ) \ {x = y}- Moreover the 
integrals 

G%x,y)dy 



jx]0,5[ 

are uniformly bounded for e g]0, 1[ and x in a compact subset of ljx]0,6[, 
and converge to 1 as e tends to 0. 

Proof By (IHTB . (lOHl) and (f6li)]l . it follows that 

GUx,y) = \detT'{z)\G'{T{z)x,T{z)y)) 

yr2(,(i_,))('^+2)/2 



detT'(z) 



(|r(z)(x-y)|2 + 4ex„y„)(«+«+2)/2 
is nonnegative. Thus 

(jya+2^ia+2)/2 



(6.37) GUx,y) + \d,GUx,y)\< 



(|a;-y|2 + ex„y„)("+"+2)2- 



In particular, G^{x,y) = Gy{x,y) is also dominated by the same bound, 
implying that G^{x,y) ^ as e ^ when x ^ y, uniformly on compact 
subsets. 

Integrating in the tangential variables first implies that 

The integral over — yn| ^ 2;„/2 tends to zero, uniformly for x„ in a 
compact subset of ]0, 5]. The integral over the remaining interval {jx^— yn| ^ 
x„/2 is bounded by 

^ (j^a+2^{a+2)/2 C'£(a+2)/2 

dyn= I 77-^^^— uTZ^dy 



^ (|x„-y„|2 + £x2)(-+3)2 7_i (|y„|2 + e){«+3)2 



5 uynr + 

- ' (|.„|2+'i)('^+3)2^^" < 



Moreover, by (j6.37p there holds 

C7|x-y|y;j+2e('^+2)/2 



\GUx,y)-Glix,y)\ 



< 



(|x-y|2 + exnyn)(''+"+^)2 
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This implies that 



Jim/ JGl{x,y) -G'y{x,y)\dy = 0. 



ix]0,5[ 

By Lemma [63] (see |GS] ). for x G there holds 



lim, / Gl{x,y)dy = 1, 

Ja;x]0,<5[ 

implying that 



lim / G'^{x,y)dy 
Jojx]O.S\ 



This finishes the proof of the lemma. □ 

We denote by E, etc the operator with kernel E, etc. In particular, 
since the kernels are smooth for ?/„ > and bounded on x M", the 
operators E*", and K'^ are defined from the space of integrable functions 
with support in w x [0,6] to the space of C°°(M_,_). The lemma above implies 
that 

(6.39) £E= = + 

Moreover, the estimates (j6.20p and (j6.35p imply that the operators E"^ [resp. 
K^] are uniformly bounded from LP{ujx]0,6[) to L^^^(R_^_) and converge in 
the strong topology to the operators E [resp. K] defined by the kernels E 
[resp. K]. 

Proposition 6.8. The operator E maps LP{ojx]0,5[) toWlf^^iWl). For all 
relatively compact open set lui C to, all S' < S and for all x £ C^{uj x [0, 6[) 
such that X = 1 on a neighborhood o/uJi x [0,5] , there is a constant C such 
that for all p G [2, +oo[ and all g £ supported in ujix]0, S'[, there holds 

(6.40) [[xnXdjdk^g[[LP{Ri) + \\xdj'E,g[[LP(R'j^) < pC[[Eg[[LP^^ny 

Moreover, for all p > n and g £ supported in uJi x [0, 6'] and all 
fi < 1 — n/p, there holds 



(6.41) CEg-g = KgG Cf 



loc\ 



Proof. By Propositions 16 . 5^ for g £ LP{ujx]0, 6[), E^g and K^(7 are uniformly 
bounded in L^(M_|_) and converge in Lf^^ to Eg and K.g respectively. 
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Suppose that g is continuous with compact support in wxJO, 5[. Then, by 
Lemma [6. 71 G'^g converges to 5 as e tends to 0. This shows that the identity 
(j6.4ip is satisfied, in the sense of distributions, for g ^ with compact 
support in a;x]0,5[. By density, this identity extends to g £ L'^{u)x]0, 6[). 
In particular. Eg € Lf^^{Wl) and CBg G Lf^^{M^). Using the regularity 
properties of L proved in [BCj . this implies that E maps L^{ujx]0,6[) into 
the space H^j^^J^u; x [0, 5[) of functions u G H^^^{uj x [0,6[) such that x„n G 

In particular, this proves that for xi ^-nd X2 in Co°(M4,), the operators 
XiXndjdkEx2 and xidjEx2 are bounded from L^(M" ) to L^(M" ). Together 
with the estimates of Proposition 16.51 and with Proposition 16.31 this implies 
that E satisfies the properties listed in the proposition and the estimates 

The identity (j6.4ip holds for g €z thus for g G L°°. The properties of 
K follow from Proposition 16 . 21 together with the estimates of Lemma [6.61 □ 



6.3 Proof of Theorem 16.11 

By assumption, cj) and g = Ccj) have compact support in x [0, 6[. Let x 
and X in Cq^{uj x [0,S[) such that x = 1 on a neighborhood of the support 
of (j) and X = 1 on a neighborhood of the support of x- 

By Propositions 16.51 and 16.21 x^9 ^ C'C(M^). Let V = X^9i which 
satisfies the estimates ()6.40p . Moreover, 

C^P = xCBg + [£, x]^g = g + X^9 + x]x^9- 

In particular, 

(6.42) C{<P-iP)=xKg+[C,x]x'^9. 

By Proposition ESI G Cf(l+) and since xEg G Cf(M_,_), the com- 

mutator [£, xlxEg also belongs to C'^(M^). Thus the right hand side h of 
([02]) is in C^'{Wl) 

We are now in position to apply the results of [GS] (or [BCM| ). which 
imply that 

<^-vecr„:\M:), x„(,^-v.)GCf,f(i:). 

Therefore, 

\\xndjdk{(l) - ip)\\L'^{Ri)+\\dj{4> - V')IIl°°(R!{:) < 

(6.43) / N 
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with h the right-hand side of (j6.42p . 

Together with the estimates of Proposition l6.8l for the LP norm of Xndjdkip 
and 9jV') this imphes the theorem. 

7 Proof of global existence, uniqueness and regu- 
larity theorem. 

Regularity of weak solutions. Suppose that {v,oj) is a weak solution 
in the sense of Definition EH Then v G L~(0, T; ^^^'^(17)) for all p < 
+00 from the elliptic regularity theorem. Therefore, using Di-Perna and 
Lions renormalization technics for transport equation p.3p . implies that 
uj G C{0,T;U{n)), hence that v G C(0, T; for all r < +00. The 

other properties of v also follow from Theorem 12.31 

Existence of a global weak solution. We construct a global weak solu- 
tion as the inviscid limit of solutions to a Navier-Stokes system with artificial 
viscosity and boundary conditions. 

Global existence for a viscous vorticity- stream function formulation. We 
define b^ = b + e and consider the following system 

' dtib^uj^) + bv' ■ Vuj^ - ediv(6^Va;") = in 17, 
6''a;^|t=o = h^'ojQ in Vt, u}^\dn = 0, 

C^-^) ' - div(^V^") = b^uj^ in ^^\an = 0, 

b 

The existence of a global weak solution [v"^ ^uj"^) for (Jq G -L^(O) follows stan- 
dard techniques for equations of Navier-Stokes type since ¥ is a strictly 
positive function, see jLQTj . From the a-priori estimates, we get the follow- 
ing uniform bounds with respect to e 

Vb'uj' G L~(0,r;L2(17)), ^/b'Vuj' G L^{Q,T;L^{Vt)). 

Moreover Vbv^ G L°°(0, T; L2(17)) uniformly 

The viscous limit e — > 0. Let us now assume that G L°°(r2). Since 
oj^ G L~(0,r;L2(Q)) n L2(0, T; i?o^(17)), using the elliptic regularity for 
degenerate elliptic equation, we get dtoj"^ G L°°(0, T; iJ~^(r2)) and thus G 
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C{[0,T);L^in)). We multiply the viscous equation by \Tr{uj'')\p-^Tr{uj') 
where Tji{u>^) = max(min(a'^, i?), —R) for R> 0, and we get 

- I h'\TR{Lo'){t)\P + e{p - I) f ds I b'\VTR{u:')\M'"^ dx 
P Jn Jo Jn 

= - [ b'\TR{u;oWdx. 
P Jn 

Therefore 

\\b'/PTRiu;)it)U,^n) < \\b'/'^TR{ 
Letting R go to +cxd, this gives 

\m'/Pu;'\\L.in) < ll(6')'/^^o||Lp(n). 
and since ujo G and the estimate is uniform 

ll'^'^llL°°{n) < ll'^o||L°°(f7)- 

Using the estimate obtained in the main estimate part we also get 

Wv^^h^n) < C(lk^llL-(n) + \\Vbv^\\L2(^n))- 

Next we note that the set {6^tJ^,e €]0, 1]} is relatively compact in 
C70([0,cx));L2 (O)) and in C^{[0,oo); L^^{Q)) as in [ LOT] . Thus {b^uj'} is 
a relatively compact set in Lf^^{[0,oo); H~^{Q)). Therefore {Vbv^} is rela- 
tively compact in Lf^^{[0, oo); L'^{Q)). This allows to pass to the limit in the 
viscous formulation and to get the global existence of a weak solution of the 
inviscid stream- vorticity formulation. 

Uniqueness of weak solutions. 

Proof. Let vi and V2 be two solutions of (jl.ip . Then v = vi — V2 satisfies 
dtv + V2 • + Vp = —V • Vvi. 

Therefore, 

(7.2) j^\\Vbv\\%<2l^b\v\'\Vvi\ 
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If Vfi G L°°, this clearly implies that v = 0. But in general, this L°° 
estimate is not available, but following [YoJ, sharp estimates can be 
substituted. By Theorem 12.31 

In the right hand side of (j7.2p . we use Young's inequality with 

\Vvi\eLP, b'/p'\v\'/p' eLP', b'-'/p'\v\^-^/p' = \bv\'/PeL^, 

implying that y{t) := ||\/6w(t)||^2 satisfies for all p > 3: 

(7.4) dty{t)<pCM^/P{y{t)}'/p' 
with M = \\Vbv\\L°° + II V^'t^2||L°°- Optimizing in p yields 

(7.5) a.y(t, < eCy(t)^^^^^^y 

Since y(0) = 0, this implies that y = 0. □ 
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